/* Copyright  2009-2010
   @Author
   Richard Changde Yin            e-mail yinchangde@hotmail.com

   This program is free software; you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation; version 2 of the License.

   This program is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the Free Software
   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA */

#ifndef _decimal_h
#define _decimal_h

#if defined(__cplusplus)
extern "C" {
#endif
#define mi_sint1korr(A) ((int8_t)(*A))
#define mi_sint2korr(A) ((int16_t) (((int16_t) (((const unsigned char*) (A))[1])) +\
                                  ((int16_t) ((int16_t) ((const char*) (A))[0]) << 8)))
#define mi_sint3korr(A) ((int32_t) (((((const unsigned char*) (A))[0]) & 128) ? \
                                  (((uint32_t) 255L << 24) | \
                                   (((uint32_t) ((const unsigned char*) (A))[0]) << 16) |\
                                   (((uint32_t) ((const unsigned char*) (A))[1]) << 8) | \
                                   ((uint32_t) ((const unsigned char*) (A))[2])) : \
                                  (((uint32_t) ((const unsigned char*) (A))[0]) << 16) |\
                                  (((uint32_t) ((const unsigned char*) (A))[1]) << 8) | \
                                  ((uint32_t) ((const unsigned char*) (A))[2])))
#define mi_sint4korr(A) ((int32_t) (((int32_t) (((const unsigned char*) (A))[3])) +\
                                  ((int32_t) (((const unsigned char*) (A))[2]) << 8) +\
                                  ((int32_t) (((const unsigned char*) (A))[1]) << 16) +\
                                  ((int32_t) ((int16_t) ((const char*) (A))[0]) << 24)))
#define mi_sint8korr(A) ((int64_t) mi_uint8korr(A))
#define mi_uint2korr(A) ((uint16_t) (((uint16_t) (((const unsigned char*) (A))[1])) +\
                                   ((uint16_t) (((const unsigned char*) (A))[0]) << 8)))
#define mi_uint3korr(A) ((uint32_t) (((uint32_t) (((const unsigned char*) (A))[2])) +\
                                   (((uint32_t) (((const unsigned char*) (A))[1])) << 8) +\
                                   (((uint32_t) (((const unsigned char*) (A))[0])) << 16)))
#define mi_uint4korr(A) ((uint32_t) (((uint32_t) (((const unsigned char*) (A))[3])) +\
                                   (((uint32_t) (((const unsigned char*) (A))[2])) << 8) +\
                                   (((uint32_t) (((const unsigned char*) (A))[1])) << 16) +\
                                   (((uint32_t) (((const unsigned char*) (A))[0])) << 24)))
#define mi_uint5korr(A) ((uint64_t)(((uint32_t) (((const unsigned char*) (A))[4])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[3])) << 8) +\
                                    (((uint32_t) (((const unsigned char*) (A))[2])) << 16) +\
                                    (((uint32_t) (((const unsigned char*) (A))[1])) << 24)) +\
                                    (((uint64_t) (((const unsigned char*) (A))[0])) << 32))
#define mi_uint6korr(A) ((uint64_t)(((uint32_t) (((const unsigned char*) (A))[5])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[4])) << 8) +\
                                    (((uint32_t) (((const unsigned char*) (A))[3])) << 16) +\
                                    (((uint32_t) (((const unsigned char*) (A))[2])) << 24)) +\
                        (((uint64_t) (((uint32_t) (((const unsigned char*) (A))[1])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[0]) << 8)))) <<\
                                     32))
#define mi_uint7korr(A) ((uint64_t)(((uint32_t) (((const unsigned char*) (A))[6])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[5])) << 8) +\
                                    (((uint32_t) (((const unsigned char*) (A))[4])) << 16) +\
                                    (((uint32_t) (((const unsigned char*) (A))[3])) << 24)) +\
                        (((uint64_t) (((uint32_t) (((const unsigned char*) (A))[2])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[1])) << 8) +\
                                    (((uint32_t) (((const unsigned char*) (A))[0])) << 16))) <<\
                                     32))
#define mi_uint8korr(A) ((uint64_t)(((uint32_t) (((const unsigned char*) (A))[7])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[6])) << 8) +\
                                    (((uint32_t) (((const unsigned char*) (A))[5])) << 16) +\
                                    (((uint32_t) (((const unsigned char*) (A))[4])) << 24)) +\
                        (((uint64_t) (((uint32_t) (((const unsigned char*) (A))[3])) +\
                                    (((uint32_t) (((const unsigned char*) (A))[2])) << 8) +\
                                    (((uint32_t) (((const unsigned char*) (A))[1])) << 16) +\
                                    (((uint32_t) (((const unsigned char*) (A))[0])) << 24))) <<\
                                    32))

/* This one is for uniformity */
#define mi_int1store(T,A) *((unsigned char*)(T))= (unsigned char) (A)

#define mi_int2store(T,A)   { uint32_t def_temp= (uint32_t) (A) ;\
                              ((unsigned char*) (T))[1]= (unsigned char) (def_temp);\
                              ((unsigned char*) (T))[0]= (unsigned char) (def_temp >> 8); }
#define mi_int3store(T,A)   { /*lint -save -e734 */\
                              uint32_t def_temp= (uint32_t) (A);\
                              ((unsigned char*) (T))[2]= (unsigned char) (def_temp);\
                              ((unsigned char*) (T))[1]= (unsigned char) (def_temp >> 8);\
                              ((unsigned char*) (T))[0]= (unsigned char) (def_temp >> 16);\
                              /*lint -restore */}
#define mi_int4store(T,A)   { uint32_t def_temp= (uint32_t) (A);\
                              ((unsigned char*) (T))[3]= (unsigned char) (def_temp);\
                              ((unsigned char*) (T))[2]= (unsigned char) (def_temp >> 8);\
                              ((unsigned char*) (T))[1]= (unsigned char) (def_temp >> 16);\
                              ((unsigned char*) (T))[0]= (unsigned char) (def_temp >> 24); }
#define mi_int5store(T,A)   { uint32_t def_temp= (uint32_t) (A),\
                              def_temp2= (uint32_t) ((A) >> 32);\
                              ((unsigned char*) (T))[4]= (unsigned char) (def_temp);\
                              ((unsigned char*) (T))[3]= (unsigned char) (def_temp >> 8);\
                              ((unsigned char*) (T))[2]= (unsigned char) (def_temp >> 16);\
                              ((unsigned char*) (T))[1]= (unsigned char) (def_temp >> 24);\
                              ((unsigned char*) (T))[0]= (unsigned char) (def_temp2); }
#define mi_int6store(T,A)   { uint32_t def_temp= (uint32_t) (A),\
                              def_temp2= (uint32_t) ((A) >> 32);\
                              ((unsigned char*) (T))[5]= (unsigned char) (def_temp);\
                              ((unsigned char*) (T))[4]= (unsigned char) (def_temp >> 8);\
                              ((unsigned char*) (T))[3]= (unsigned char) (def_temp >> 16);\
                              ((unsigned char*) (T))[2]= (unsigned char) (def_temp >> 24);\
                              ((unsigned char*) (T))[1]= (unsigned char) (def_temp2);\
                              ((unsigned char*) (T))[0]= (unsigned char) (def_temp2 >> 8); }
#define mi_int7store(T,A)   { uint32_t def_temp= (uint32_t) (A),\
                              def_temp2= (uint32_t) ((A) >> 32);\
                              ((unsigned char*) (T))[6]= (unsigned char) (def_temp);\
                              ((unsigned char*) (T))[5]= (unsigned char) (def_temp >> 8);\
                              ((unsigned char*) (T))[4]= (unsigned char) (def_temp >> 16);\
                              ((unsigned char*) (T))[3]= (unsigned char) (def_temp >> 24);\
                              ((unsigned char*) (T))[2]= (unsigned char) (def_temp2);\
                              ((unsigned char*) (T))[1]= (unsigned char) (def_temp2 >> 8);\
                              ((unsigned char*) (T))[0]= (unsigned char) (def_temp2 >> 16); }
#define mi_int8store(T,A)   { uint32_t def_temp3= (uint32_t) (A),\
                              def_temp4= (uint32_t) ((A) >> 32);\
                              mi_int4store((unsigned char*) (T) + 0, def_temp4);\
                              mi_int4store((unsigned char*) (T) + 4, def_temp3); }

#ifdef WORDS_BIGENDIAN

#define mi_float4store(T,A) { ((unsigned char*) (T))[0]= ((const unsigned char*) &A)[0];\
                              ((unsigned char*) (T))[1]= ((const unsigned char*) &A)[1];\
                              ((unsigned char*) (T))[2]= ((const unsigned char*) &A)[2];\
                              ((unsigned char*) (T))[3]= ((const unsigned char*) &A)[3]; }

#define mi_float4get(V,M)   { float def_temp;\
                              ((unsigned char*) &def_temp)[0]= ((const unsigned char*) (M))[0];\
                              ((unsigned char*) &def_temp)[1]= ((const unsigned char*) (M))[1];\
                              ((unsigned char*) &def_temp)[2]= ((const unsigned char*) (M))[2];\
                              ((unsigned char*) &def_temp)[3]= ((const unsigned char*) (M))[3];\
                              (V)= def_temp; }

#define mi_float8store(T,V) { ((unsigned char*) (T))[0]= ((unsigned char*) &V)[0];\
                              ((unsigned char*) (T))[1]= ((unsigned char*) &V)[1];\
                              ((unsigned char*) (T))[2]= ((unsigned char*) &V)[2];\
                              ((unsigned char*) (T))[3]= ((unsigned char*) &V)[3];\
                              ((unsigned char*) (T))[4]= ((unsigned char*) &V)[4];\
                              ((unsigned char*) (T))[5]= ((unsigned char*) &V)[5];\
                              ((unsigned char*) (T))[6]= ((unsigned char*) &V)[6];\
                              ((unsigned char*) (T))[7]= ((unsigned char*) &V)[7]; }

#define mi_float8get(V,M)   { double def_temp;\
                              ((unsigned char*) &def_temp)[0]= ((unsigned char*) (M))[0];\
                              ((unsigned char*) &def_temp)[1]= ((unsigned char*) (M))[1];\
                              ((unsigned char*) &def_temp)[2]= ((unsigned char*) (M))[2];\
                              ((unsigned char*) &def_temp)[3]= ((unsigned char*) (M))[3];\
                              ((unsigned char*) &def_temp)[4]= ((unsigned char*) (M))[4];\
                              ((unsigned char*) &def_temp)[5]= ((unsigned char*) (M))[5];\
                              ((unsigned char*) &def_temp)[6]= ((unsigned char*) (M))[6];\
                              ((unsigned char*) &def_temp)[7]= ((unsigned char*) (M))[7]; \
                              (V)= def_temp; }
#else

#define mi_float4store(T,A) { ((unsigned char*) (T))[0]= ((unsigned char*) &A)[3];\
                              ((unsigned char*) (T))[1]= ((unsigned char*) &A)[2];\
                              ((unsigned char*) (T))[2]= ((unsigned char*) &A)[1];\
                              ((unsigned char*) (T))[3]= ((unsigned char*) &A)[0]; }

#define mi_float4get(V,M)   { float def_temp;\
                              ((unsigned char*) &def_temp)[0]= ((unsigned char*) (M))[3];\
                              ((unsigned char*) &def_temp)[1]= ((unsigned char*) (M))[2];\
                              ((unsigned char*) &def_temp)[2]= ((unsigned char*) (M))[1];\
                              ((unsigned char*) &def_temp)[3]= ((unsigned char*) (M))[0];\
                              (V)= def_temp; }

#if defined(__FLOAT_WORD_ORDER) && (__FLOAT_WORD_ORDER == __BIG_ENDIAN)
#define mi_float8store(T,V) { ((unsigned char*) (T))[0]= ((unsigned char*) &V)[3];\
                              ((unsigned char*) (T))[1]= ((unsigned char*) &V)[2];\
                              ((unsigned char*) (T))[2]= ((unsigned char*) &V)[1];\
                              ((unsigned char*) (T))[3]= ((unsigned char*) &V)[0];\
                              ((unsigned char*) (T))[4]= ((unsigned char*) &V)[7];\
                              ((unsigned char*) (T))[5]= ((unsigned char*) &V)[6];\
                              ((unsigned char*) (T))[6]= ((unsigned char*) &V)[5];\
                              ((unsigned char*) (T))[7]= ((unsigned char*) &V)[4];}

#define mi_float8get(V,M)   { double def_temp;\
                              ((unsigned char*) &def_temp)[0]= ((unsigned char*) (M))[3];\
                              ((unsigned char*) &def_temp)[1]= ((unsigned char*) (M))[2];\
                              ((unsigned char*) &def_temp)[2]= ((unsigned char*) (M))[1];\
                              ((unsigned char*) &def_temp)[3]= ((unsigned char*) (M))[0];\
                              ((unsigned char*) &def_temp)[4]= ((unsigned char*) (M))[7];\
                              ((unsigned char*) &def_temp)[5]= ((unsigned char*) (M))[6];\
                              ((unsigned char*) &def_temp)[6]= ((unsigned char*) (M))[5];\
                              ((unsigned char*) &def_temp)[7]= ((unsigned char*) (M))[4];\
                              (V)= def_temp; }

#else
#define mi_float8store(T,V) { ((unsigned char*) (T))[0]= ((unsigned char*) &V)[7];\
                              ((unsigned char*) (T))[1]= ((unsigned char*) &V)[6];\
                              ((unsigned char*) (T))[2]= ((unsigned char*) &V)[5];\
                              ((unsigned char*) (T))[3]= ((unsigned char*) &V)[4];\
                              ((unsigned char*) (T))[4]= ((unsigned char*) &V)[3];\
                              ((unsigned char*) (T))[5]= ((unsigned char*) &V)[2];\
                              ((unsigned char*) (T))[6]= ((unsigned char*) &V)[1];\
                              ((unsigned char*) (T))[7]= ((unsigned char*) &V)[0];}

#define mi_float8get(V,M)   { double def_temp;\
                              ((unsigned char*) &def_temp)[0]= ((unsigned char*) (M))[7];\
                              ((unsigned char*) &def_temp)[1]= ((unsigned char*) (M))[6];\
                              ((unsigned char*) &def_temp)[2]= ((unsigned char*) (M))[5];\
                              ((unsigned char*) &def_temp)[3]= ((unsigned char*) (M))[4];\
                              ((unsigned char*) &def_temp)[4]= ((unsigned char*) (M))[3];\
                              ((unsigned char*) &def_temp)[5]= ((unsigned char*) (M))[2];\
                              ((unsigned char*) &def_temp)[6]= ((unsigned char*) (M))[1];\
                              ((unsigned char*) &def_temp)[7]= ((unsigned char*) (M))[0];\
                              (V)= def_temp; }
#endif /* __FLOAT_WORD_ORDER */
#endif /* WORDS_BIGENDIAN */

typedef enum
{TRUNCATE=0,
 HALF_EVEN,
 HALF_UP,
 CEILING,
 FLOOR}
  decimal_round_mode;
typedef int32_t decimal_digit_t;

typedef struct st_decimal_t {
  int    intg, frac, len;
  bool sign;
  decimal_digit_t *buf;
} decimal_t;

int internal_str2dec(char *from, decimal_t *to, char **end,
                     bool fixed);
int decimal2string(decimal_t *from, char *to, int *to_len,
                   int fixed_precision, int fixed_decimals,
                   char filler);
int decimal2uint64_t(decimal_t *from, uint64_t *to);

int uint64_t2decimal(uint64_t from, decimal_t *to);

int decimal2int64_t(decimal_t *from, int64_t *to);

int int64_t2decimal(int64_t from, decimal_t *to);

int decimal2double(decimal_t *from, double *to);

int double2decimal(double from, decimal_t *to);

int decimal_actual_fraction(decimal_t *from);

int decimal2bin(decimal_t *from, unsigned char *to, int precision, int scale);

int bin2decimal(const unsigned char *from, decimal_t *to, int precision, int scale);

int decimal_size(int precision, int scale);
int decimal_bin_size(int precision, int scale);
int decimal_result_size(decimal_t *from1, decimal_t *from2, char op,
                        int param);

int decimal_intg(decimal_t *from);
int decimal_add(decimal_t *from1, decimal_t *from2, decimal_t *to);
int decimal_sub(decimal_t *from1, decimal_t *from2, decimal_t *to);
int decimal_cmp(decimal_t *from1, decimal_t *from2);
int decimal_mul(decimal_t *from1, decimal_t *from2, decimal_t *to);
int decimal_div(decimal_t *from1, decimal_t *from2, decimal_t *to,
                int scale_incr);
int decimal_mod(decimal_t *from1, decimal_t *from2, decimal_t *to);
int decimal_round(decimal_t *from, decimal_t *to, int new_scale,
                  decimal_round_mode mode);
int decimal_is_zero(decimal_t *from);
void max_decimal(int precision, int frac, decimal_t *to);

#define string2decimal(A,B,C) internal_str2dec((A), (B), (C), 0)
#define string2decimal_fixed(A,B,C) internal_str2dec((A), (B), (C), 1)

/* set a decimal_t to zero */

#define decimal_make_zero(dec)        do {                  \
                                        (dec)->buf[0]=0;    \
                                        (dec)->intg=1;      \
                                        (dec)->frac=0;      \
                                        (dec)->sign=0;      \
                                      } while(0)

/*
  returns the length of the buffer to hold string representation
  of the decimal (including decimal dot, possible sign and \0)
*/

#define decimal_string_size(dec) (((dec)->intg ? (dec)->intg : 1) + \
                                 (dec)->frac + ((dec)->frac > 0) + 2)

/* negate a decimal */
#define decimal_neg(dec) do { (dec)->sign^=1; } while(0)

/*
  conventions:

    decimal_smth() == 0     -- everything's ok
    decimal_smth() <= 1     -- result is usable, but precision loss is possible
    decimal_smth() <= 2     -- result can be unusable, most significant digits
                               could've been lost
    decimal_smth() >  2     -- no result was generated
*/

#define E_DEC_OK                0
#define E_DEC_TRUNCATED         1
#define E_DEC_OVERFLOW          2
#define E_DEC_DIV_ZERO          4
#define E_DEC_BAD_NUM           8
#define E_DEC_OOM              16

#define E_DEC_ERROR            31
#define E_DEC_FATAL_ERROR      30

#if defined(__cplusplus)
}
#endif

#endif

